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Understanding the dynamics of ions in nanopores is essential for applications ranging from single- 
molecule detection to DNA sequencing. We show both analytically and by means of molecular 
dynamics simulations that under specific conditions ion-ion interactions in nanopores lead to the 
phenomenon of ionic Coulomb blockade, namely the build-up of ions inside a nanopore with specific 
capacitance impeding the fiow of additional ions due to Coulomb repulsion. This is the counterpart 
of electronic Coulomb blockade observed in mesoscopic systems. We discuss the analogies and 
differences with the electronic case as well as experimental situations in which this phenomenon 
could be detected. 



Recently, there has been a surge of interest in 
nanopores due to their potential in several technologi- 
cal applications, the most notable being DNA sequencing 
and detection [TH7]. In addition, nanopores offer unprece- 
dented opportunities to study several fundamental physi- 
cal processes related to ionic transport in confined geome- 
tries. For example, due to the tightly bound hydration 
layers surrounding an ion, it may be possible to observe 
steps in the ionic conductance due to the shedding of 
water molecules in the individual layers as the ion-water 
"quasi-particle" passes through the pore opening |H| 
This can be viewed as the classical analog of the electron 
conductance quantization in nanoscopic/mesoscopic sys- 
tems (see, e.g., Ref. [10, 11 ). 

In this letter, we discuss another fundamental physical 
process which may occur in nanopores due to screened 
ion-ion interactions. This many-body effect originates 
when ions flow under the effect of an electric (or pressure) 
field and build up in nanopores of specific capacitances, 
so that the flow of additional ions is hindered. To be 
specific. Fig. [T| shows a typical pore geometry we have in 
mind with an opening wider than the other so that ions of 
one polarity can easily accumulate inside the pore, while 
ions of opposite polarity will accumulate outside the nar- 
row opening of the pore. The effect we consider would 
then occur when the kinetic energy of an ion in solution, 
given approximately by the sum of the thermal energy 
and M(/x£')^/2 (where n is the ion mobility, M its mass 
and E the electric field) is smaller than the capacitive 
energy /2C, with Q the average number of ions in the 
pore, and C the capacitance of the pore. 

This effect is similar to the Coulomb blockade phe- 
nomenon experienced by electrons transporting through 
a quantum dot weakly coupled to two electrodes - in the 
sense that the tunneling resistance between the electrodes 
and the dot is much larger than the quantum of resis- 
tance |12| . Like the quantum case, we also expect ionic 
Coulomb blockade to occur when the nanopore is "weakly 
coupled" to at least one ion reservoir. In this case, this 
means that the average time it takes an ion to "cross" the 
weak contact link is much longer than the average time it 



takes the same ion to propagate unhindered by the pore. 
In other words, the "contact" ionic resistance of the pore 
with at least one reservoir has to be much larger than 
the resistance offered to ion flow without the pore. From 
previous work [8, 9J we know that such a situation may 
be realized when the pore opening is in the "quantum" 
regime, namely when the pore radius is small enough that 
an ion has to shed part of its hydration layers to cross 
the pore aperture. 

The analogy with the electron case therefore seems 
complete safe for a notable exception. In the ionic case, 
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Figure 1. (Color online) Right panel: a snapshot of a molec- 
ular dynamics simulation of IM KCl solution subject to an 
electric field and translocating through a cone-shaped Si3N4 
pore (indicated by the yellow and blue atoms) at a pore slope 
angle of 45° with a bottom opening of 7 A radius and thickness 
of 25 A. Cl~ ions and K"^ ions are colored aqua and brown, 
respectively. The symbols 74 and 7;, indicate the rates of ion 
transfer at the bottom and top openings, respectively. The 
field used to generate this figure is 5.0 kcal/(molAe) which 
allows an easier visualization of the ions build-up. Left panel: 
the net charge density corresponding to the right panel config- 
uration. It is clear that in this case there is an accumulation of 
Cl~ ions inside the pore while the K""" ions are located mostly 
outside the pore. 
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the Pauli exclusion principle need not be satisfied: we can 
fit as many ions inside the pore as its capacitance allows 
without limitations from Fermi statistics. We thus expect 
ionic Coulomb blockade to occur as a function of the ionic 
concentration - which could be thought of as controlled 
by a "gate voltage" - but not as a function of bias, namely 
we expect the current to show a strongly non-linear be- 
havior as a function of molarity - within the limits of ion 
precipitation - at fixed voltage, and essentially linear be- 
havior with bias - for values below electrolysis - at fixed 
molarity [13) . In addition, nanopores as those shown in 
Fig. [ijcan be easily fabricated [T]. We thus expect our 
predictions to be readily accessible experimentally. 

Note also that the phenomenon we consider here is 
fundamentally different from other nonlinear ionic trans- 
port effects that have been discussed previously in the 
literature [51 [T^HTT] . For instance, ions in solution can 
form an "ionic atmosphere", namely a region around a 
given ion in which ions of opposite charge are attracted 
electrostatically [14,. However, the ionic atmosphere is a 
dynamic phenomenon with typical lifetimes on the order 
of 10^^ s and it can be easily destroyed at field strengths 
of 10"* V/cm ^14]. Therefore, for the fields and molarities 
we consider in this work the ionic atmosphere is not rel- 
evant |14| . Another deviation from Ohmic behavior has 
been observed experimentally in Ref. |15j. In this work, 
the addition of a small amount of divalent cations to an 
ionic solution leads to current oscillations and negative- 
incremental resistance. The phenomenon has been ex- 
plained as due to the transient formation and breakup of 
nano-precipitates |15| . The nano-precipitates temporar- 
ily block the ionic current through the pore thus causing 
current oscillations. Again, this effect does not pertain 
to the present work. 

Let us then start by presenting first a simple analytical 
model that captures the main physics of ionic Coulomb 
blockade. We will later corroborate these results with 
all-atom classical molecular dynamics simulations. Like 
the electronic Coulomb blockade case (see, e.g., [TS]), we 
can model the dynamics of ionic conduction using a rate 
equation approach by calculating the transition proba- 
bilities for the pore to accommodate a certain amount of 
charges in addition to the background charge. Referring 
again to Fig. [TJwe consider the rates of conduction in and 
out of the wide opening of the pore ("top opening") and 
the neck of the pore ("bottom opening") , as indicated by 
7t and 7b, respectively. For clarity, we refer to the ions 
with the charge sign that favors accumulation inside the 
pore. In the case of Fig. [T] these are Cl~ ions, as also 
demonstrated by the density profile on the left panel of 
the same figure. The opposite-charge ions accumulate on 
the outskirts of the pore entrance (see Fig. [ij and the 
current blockade occurs when ions from the reservoir ad- 
jacent to the neck of the pore attempt to enter it. The 
considerations we make below would then be similar, but 
clearly with different values of the parameters and the di- 



rection of ion motion. If we reversed the bias direction we 
would obtain similar results but with the opposite charge 
sign accumulation inside and outside the pore. 

We assume that the rate through the bottom opening 
is related to the top one by 



lb = ajt 



< a < 1 , 



with 



7t = At^iUoE , 



(1) 



(2) 



where At is the top area of the pore, /i is the ionic mo- 
bility, no is the ionic density, and E = V/d is the electric 
field, where V is the voltage and d is the length of the 
pore |19j . The parameter a takes into account the dif- 
ference in the two rates due to the partial shedding of 
the hydration layers around the ions when they cross the 
bottom neck, and the fact that the radius of the pore 
opening is smaller at the neck than at the top. 

In order to be able to solve the rate equations analyti- 
cally we assume only two possible "ionic states". The first 
state with transition probability Pq corresponds to only 
the background charge in the pore, namely to [npfipj 
ions, where rig is the background density, Qp is the vol- 
ume of the pore, and the symbol [■ • • J represents the 
fioor function which will make the product npilp an in- 
teger. The second state with probability Pi is that of 
the background charge plus one single extra ion, namely 
[noilpj -I- 1 ions in the pore. While this may seem like 
an oversimplified situation, it captures the main trends 
observed with the molecular dynamics simulations as we 
show below. The Markovian equation for the transition 
probability Pq is then 



dl^ 
dt 



(3) 



where Fo-i-i is the transition rate of going from state to 
state 1, and Fi^o is the reverse process. The equation of 
motion for Pi is easily obtained by interchanging with 
1 in Eq. [S) 

The transition rate Fo-i-i contains the capacitive bar- 
rier experienced by an ion entering the pore already oc- 
cupied by [rioilpj ions. Using the steady-state Nernst- 
Planck equation it can then be written as [5] 



7f exp(-ec/fcr), 



(4) 



where ec is a single-particle capacitive energy barrier, k 
the Boltzmann constant, and T the temperature. We can 
write the single-particle capacitive energy as jl8| 



2[nof^pJ + 1 
2C 



(5) 



where e is the elementary charge and C is the capacitance 
of the pore. Since the capacitance of the pore is linearly 
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related to the ratio of the surface areas of the top and 
bottom openings, it must vary according to the param- 
eter a we have introduced in Eq. ([ij. We then assume 
C = aCo with Cq some reasonable experimental value 
for the capacitance. The rate Ti^q is instead simply 



76, 



(6) 



because there is no capacitive barrier for this process. We 
also note that, in principle, there is a finite rate for an 
ion to move in the direction opposite to that determined 
by the electric field. This rate would contribute to both 
Fo-i-i and ri_>.o. However, this process is exponentially 
suppressed for the biases we consider here, and we can 
thus safely ignore it [50' . 

At steady state, dP/dt = 0, and from Eq. ^ - and the 
equivalent equation for Pi - we obtain the steady-state 
probabilities 
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Pi 



7texp(-ec/A:T) + 7^ 



7texp(-ec/A:T) 



To 



+ ri^o 7* exp (-ec/fcr) -I- 7b 



(7) 



(8) 



Finally, the current at steady state is the same every- 
where and we evaluate it across the neck of the pore 



(9) 



where the superscript b indicates that we retain only the 
terms corresponding to the bottom part of the transition 
rates. In the present case, this gives 



If, 



exp {~ec/kT) 
exp {—ec/kT) + a 



(10) 
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Figure 2. (Color online) Current from Eq. ( 10 1 as a function 
of the ion concentration at a fixed voltage of = 1.0 V for 
various pore neck openings as represented by the parameter 
a (larger a implies larger neck opening). The inset shows the 
current, at a concentration of 1.0 M, as a function of bias for 
a = 0.2. See text for all other parameters. 
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Figure 3. (Color online) Current as a function of KCl molarity 
for various neck radii and at a fixed bias of 1.0 V. For relevant 
experimental molarities, we observe an almost linear increase 
of the current for a radius r = 10 and a non-linear behavior 
approaching a saturation at a given molarity for smaller radii. 
The current saturation and decrease is shown explicitly in the 
left inset for a 7 A neck radius pore. We also show the increase 
in current with voltage in the other inset. 



with 7( and ec given by Eqs. ^ and ([S]), respectively. 
This equation shows indeed what we were expecting: the 
current is linearly dependent on bias (via the parame- 
ter 7t) for a fixed ion concentration, and it saturates - 
and eventually decreases - with increasing concentration 
according to the difference between the top and bottom 
openings (parameter a). This is shown in Fig. |2] where 
we have used the parameters At = 7.8 nm^^, d ~ 25 A, 
^ = 7 X 10-8 m^V-^s-i, T = 295 K, and C/a = 1.0 fF. 
The latter is a reasonable value for a nanopore of our 
size jH]. With these parameters we indeed find that the 
typical drift contribution to the kinetic energy of an ion 
in solution is some fraction of 1 meV, which (as the ther- 
mal energy) is smaller than the capacitive barrier (on the 
order of several tens of meVs) . 

We now show that this model captures the main 
physics of this phenomenon by performing all-atom 
molecular dynamics simulations using NAMD2 |[22j. The 
pores are made of 25 A thick silicon nitride material in 
the /3-phase and they have a conical shape with a slope 
ang le of 45° (see Fig. IT}. We vary the bottom opening 
from a radius r — 5 A to a radius r = 10 A. We then 
introduce a given concentration of KCl while keeping the 
temperature fixed at 295 K. By introducing an external 
constant electric field we can then probe the ionic con- 
ductance. Other details of the simulations are reported 
in Ref. EU. 

The results of these simulations are plotted in Fig. |3] 
As predicted by Eq. ( 10 1 our calculations show an almost 
linear behavior of the current as a function of bias, and 



for a fixed bias a saturation of the current as a func- 
tion of ion concentration which is more pronounced for 
pores with smaller neck radius. In the inset of Fig. [3] we 
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also show explicitly current saturation and decrease for a 
7 A neck radius pore, which has to be compared with the 
results reported in Fig. |2] Note, however, that current 
saturation and decrease occur at molarities well above 
the ionic precipitation limit of KCl of about 3.5 M. We 
thus expect these two features not to be directly visible 
for this configuration. 

Finally, we discuss conditions under which the phe- 
nomenon of ionic Coulomb blockade may be better re- 
solved. It is definitely true that nanopores of a V-shape 
as that shown in Fig.[l]can be easily fabricated [T]. How- 
ever, the amount of surface charges on the internal walls 
of the pore is not so straightforward to control. Never- 
theless, surface charges are beneficial to the observation 
of this effect. Indeed, if a certain amount of, say, positive 
charges are present at the surface of the pore walls, more 
Cl^ ions would accumulate inside the pore even with- 
out the need of an electric field. Everything else being 
equal, the predicted non-linearities as a function of mo- 
larity would then occur at lower concentrations. Better 
yet would be if these surface charges could be modified 
electrically. This has been recently achieved experimen- 
tally by placing electrodes inside nanopores of a shape 
similar to the one considered in this paper |231 - l25| . The 
predicted phenomenon is therefore within reach of these 
experiments. Note also that the same phenomenon would 
occur in the case of a straight- wall (e.g., cylindrical) pore 
provided the openings are small enough that the hydra- 
tion layers need to be partly shed. In this case, using 
molecular dynamics simulations we have also observed 
ionic current blockade as a function of molarity (not re- 
ported in this paper). Here, however, we expect the pore 
wall charges to play a much more important role than in 
the case of V-shaped pores. In this respect, experiments 
with more well-controlled synthetic pores, e.g., nanotube 
pores PSI may be particularly suitable to verify this phe- 
nomenon. We leave all these studies for future work. 

In conclusion, we have predicted that ionic Coulomb 
blockade effects should be observable in nanopores of spe- 
cific geometries, when ions accumulate in the pore thus 
impeding the fiow of like-charge ions. We have corrobo- 
rated these predictions with all-atom molecular dynamics 
simulations. The ensuing non-linearities as a function of 
ion concentration are within reach of experimental veri- 
fication. This phenomenon parallels the one observed in 
electronic transport across quantum dots weakly coupled 
to two electrodes. This effect may play a significant role 
in various applications of nanopores presently pursued, 
and is of importance to understand ionic conduction in 
confined geometries. 
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